Mathematical Biosciences and Engineering
● American Institute of Mathematical Sciences (AIMS)
Preprints posted in the last 7 days, ranked by how well they match Mathematical Biosciences and Engineering's content profile, based on 23 papers previously published here. The average preprint has a 0.03% match score for this journal, so anything above that is already an above-average fit.
Tuffaha, M. Z.; Dick, D. W.; Li, J.; Heffernan, J. M.
Show abstract
We develop an age- and immunity-structured COVID-19 model with setting-specific contacts, calibrated to Ontario incidence and hospital admissions before Omicron. We simulate counterfactual vaccination strategies and contact restrictions in schools, workplaces, and other community settings, evaluating projected severe-case demand against an illustrative healthcare-capacity benchmark. Vaccination combined with workplace and community contact restrictions keeps simulated demand within capacity. In contrast, relaxing restrictions in either setting produces large surges and repeated capacity threats, even when school restrictions remain. School-only changes produce smaller effects on severe disease. These findings highlight the value of vaccination and setting-specific contact restrictions that prioritize workplace and community settings for pathogens with adult-concentrated severe outcomes.
Habib, G.; Thakur, S.; Shaban, M.; Ahmad, R.; Malik, I. A.
Show abstract
The -Bonacci weighted sequence space {ell}(,) is introduced for [≥] 2, 1 [≤] < {infty}, and > 0). It is proved to be a Banach space which is isometrically isomorphic to {ell}, with an unconditional Schauder basis. The space is also reflexive and uniformly convex. Its strict inclusion chain, continuous dual, and asymptotic truncation estimate of order - / are also established. Based on this framework, Scaled 𝑚-BonacciNet is proposed as a lightweight CNN that rescales Fibonacci-Net channels using the dominant 𝑚-Bonacci root and employs depth wise separable convolutions. Tribonacci-Net achieved 96.09% accuracy and a 0.9581 macro F1-score on the 15-class Brain MRI dataset using 0.85 million parameters and 0.366 GFLOPs. Pentanacci-Net attained 95.87% accuracy on the 44-class dataset while reducing parameters and FLOPs by 47.10% and 44.48%, respectively. On the CKD ultrasound dataset, Hexanacci-Net achieved 83.75% accuracy and a 0.7941 macro F1-score, with corresponding reductions of 50.37% and 45.10%. These findings support dominant-root scaling as a structured approach to lightweight CNN design.
Sandor, D.; Nemeth, A.; Daroczy, B.; Sztanek, F.; Harangi, M.; Antal, P.
Show abstract
Personalized medicine aims to tailor treatments to individual patients, but its implementation is limited by scarce clinical resources and heterogeneous treatment effects. Electronic health records now offer longitudinal data that can support data-driven therapy selection. We focus on SGLT2 inhibitors, medications with broad metabolic and cardiovascular benefits, whose clinical impact varies significantly based on individual patient profiles. We propose a counterfactual Structural Causal Model (SCM) to estimate individual treatment effects (ITE) of SGLT2 inhibitors from observational data. The model assumes counterfactual stability and uses a Gumbel-Max sampling mechanism to capture stochastic treatment assignment while enabling differentiable inference. Optimal treatments are defined by maximizing the counterfactual value of future diagnostic outcomes. We evaluate the Counterfactual Gumbel-Max SCM using a doubly robust value estimator and show that, across diagnostic categories, it outperforms real-world prescriptions in predicted treatment value, highlighting subpopulations where SGLT2 inhibitors yield superior outcomes. By integrating causal inference with electronic health data, this framework provides individualized counterfactual reasoning about SGLT2 treatment effects and can support clinicians in identifying optimal therapies for type 2 diabetes and comorbid conditions.
Ko, Y.; Santos, C.; Huang, A. T.; Salje, H.; Cummings, D. A.; Ko, A. I.
Show abstract
Brazil experienced its worst dengue epidemic in 2024, with an increase in severe outcomes among elderly individuals, highlighting the need for vaccination strategies informed by population immunity. We reconstructed state- and age-specific immune profiles by fitting an age-structured catalytic model to dengue surveillance data from 2000-2025. This inferred immune profile was incorporated into an SIR-type transmission model to project disease burden and to evaluate the public health impact of vaccine allocation strategies. We observed substantial spatial heterogeneity in accumulated immunity, which was closely associated with projected burden and optimal allocation. Strategies prioritizing states with lower historical incidence and larger susceptible populations yielded greater benefits. Under limited vaccine supply, targeting children resulted in the greatest reductions in DALYs, whereas targeting elderly populations achieved smaller DALY reductions but larger reductions in mortality. We further showed that accelerated vaccine supply markedly increased impact, while allocation choices remained critical when availability was constrained.
Jeong, I.; Lee, T.; Park, J.-H.; Kim, Y.; Lee, H.
Show abstract
Background: Conventional discrimination metrics, such as the area under the receiver operating characteristic curve (AUROC), do not characterize the composition of outcomes among patients flagged as high risk by clinical prediction models. This limitation is particularly relevant in critical care, where competing deaths and incomplete follow-up can complicate the interpretation of prediction alerts. We developed the Balanced Risk Integrity Index (BRII), a stratum-level metric defined as the proportion of target events among target events, competing deaths, and early-censored observations in a flagged high-risk stratum. BRII is algebraically related to positive predictive value (PPV) and is reported alongside its component decomposition to characterize the outcome composition of flagged patients. Methods: We evaluated BRII using seven acute kidney injury (AKI) prediction models developed in MIMIC-IV and externally validated in MIMIC-III and eICU-CRD (116,010 patients in total). Additional simulations examined the effects of competing deaths and early censoring on BRII and conventional model performance metrics. Results: For 48-hour AKI prediction in the top 10% risk stratum, BRII ranged from 0.90 to 0.96 in MIMIC-III, compared with less than 0.26 across all models in eICU-CRD, where 61.5% of evaluated patients were censored before the prediction horizon. In the eICU-CRD top 10% strata for the same target, 37.0% to 50.2% of patients were censored early, 12.2% to 15.5% had AKI, and 3.4% to 8.1% died before AKI. BRII and PPV ordered at least one model pair differently in 34 of 48 cohort-outcome-top-k groups. Under simulated random censoring for 48-hour AKI prediction, BRII declined by up to 0.80 while binary-label AUROC changed by less than 0.02. Conclusions: BRII provided complementary information on high-risk stratum composition, yielding different model rankings from PPV and revealing substantial differences across external validation cohorts. Interpreted alongside PPV and its component decomposition, BRII offers a standardized approach to evaluating clinical prediction alerts across models and observation environments.
Westfall, M. N.; Hoiby, N.; Rasmussen, S.; Ziock, H.-J.
Show abstract
To deal with the COVID-19 pandemic, Denmark implemented an aggressive program aimed at controlling the pandemic which included both testing and contact tracing programs. The large amount of data collected in Denmark on the pandemic and its integrated nature made it possible to investigate the effectiveness of the testing and contact tracing effort. The results of that investigation and a discussion of the effectiveness and limitations of the different control efforts are reported herein, as are lessons learned concerning possible improvements for the future. The analysis is limited to the period from the start of the COVID-19 pandemic in Denmark (late Feb. 2020) into May 2021, at which time the impact of vaccination intervention was already playing a role with ~25% of the Danish population at least partially vaccinated. From the available data sets we estimate the national Infection Detection Efficiency (IDE) starting around July 1, 2020, the time at which the national infection detection program started to be implemented in a more robust form. The results rely on daily COVID-19 hospital admissions, PCR test results, and anti-body prevalence values for the Danish population. We find that for the long-term average, the combined Danish national testing and contact tracing program identified on average about 55 to 60% of the infected individuals during the period from Sept. 1, 2020, into early May 2021. However, higher infection detection efficiencies were achieved in localized areas for short time intervals. Our investigations further indicate that the Danish testing and contact tracing program, which was among the most ambitious in the world, could have had a better balance between voluntary probing on the one hand and focused contact tracing and targeted testing on the other hand. As typically implemented, it proved to be largely ineffective in controlling the pandemic. However, in at least one case when it was implemented in an enhanced and focused mode, it was able to bring a large local outbreak that dominated the country's cases under control. We have compiled a list of lessons learned regarding the use of contact tracing procedures to control potential future pandemics before vaccinations become available.
Couper, L. I.; Camponuri, S. K.; Sykes, J. E.; Li, X.; Hubbard, A. E.; Bustamante, B. L. M.; Williamson, T.; Ruberto, I.; Sondermeyer-Cooksey, G.; Karmarkar, E. N.; Yu, A. T.; Heaney, A. K.; Head, J. R.; Remais, J. V.
Show abstract
Importance: Coccidioidomycosis, also known as Valley fever, is a serious fungal disease caused by inhaling Coccidioides spores aerosolized from the soil. Despite rising public health concern, the distribution and incidence of coccidioidomycosis in the US remain poorly characterized as not all states are required to report cases, and case under-ascertainment is presumed to be high, particularly outside of the endemic region. Objective: To estimate monthly, county-level coccidioidomycosis incidence across the US while accounting for widespread under-ascertainment. Design, Setting, and Participants: We integrated multisource data, including human cases diagnosed in electronic health records, serological testing of dogs, and ecological and sociodemographic predictors, in a machine learning model estimating monthly, county-level coccidioidomycosis incidence. The model was trained on surveillance data from Arizona and California from 2013 to 2023, then applied to predict incidence in states across the contiguous US over the same time period. Main Outcomes and Measures: Incidence of coccidioidomycosis at the monthly county level. Results: The model predicted coccidioidomycosis incidence from 2013 to 2023 for 3,108 counties from all states in the contiguous US. The model identified 60 counties in non-reporting states where predicted incidence was [≥]2 cases per 100,000 population (considered moderate to high incidence). In states with mandatory reporting outside of Arizona and California, the model predicted county-level incidence rates that were 28% higher (95% bootstrap interval, 0.1%-69.8%) than those captured in national surveillance, equating to a total of 608 potentially missed cases per year (95% bootstrap interval, 340-964). Conclusions and Relevance: In this cross-sectional study, the model quantified case under-ascertainment for coccidioidomycosis relative to Arizona and California, and highlighted the potential for alternative data sources to fill in gaps in mandatory reporting. Estimating incidence while accounting for under-ascertainment and non-reporting may help efforts to improve public and physician awareness and to detect emergence in new regions.
Bou Ali, M. N.; Lassoued, N.; Ben Amor, B.; Nasr, S.; Aouididi, M. A.
Show abstract
Background. In people with type 2 diabetes (T2D), postprandial glucose weighs on glycaemic control, and carbohydrate counting is difficult. Patients can describe their meal but rarely quantify it. Objective. To assess (1) whether a small open large language model (LLM), fine-tuned on public data, estimates the carbohydrate content of a meal described in words, and (2) what this carbohydrate adds to the prediction of the postprandial glucose rise in two public cohorts. Methods. Step 1: the Qwen3-4B LLM was fine-tuned with QLoRA on Nutrition5k dishes and compared with an automatically matched USDA food composition table and with Qwen3-4B and Qwen3-8B used without training, on test dishes (internal control) and on 100 real meals of ShanghaiT2DM patients, whose reference carbohydrate was computed from the recorded quantities. Step 2: the maximum glucose rise between 0 and 120 minutes ({Delta}Gmax) was predicted by personalised gradient boosting (M2), compared with a spline linear regression and with the Chronos-2 foundation model, using only information available at meal start and person-level cross-validation, in ShanghaiT2DM (people with T2D recruited in two Shanghai hospitals; carbohydrate estimated by the LLM) and CGMacros (free-living adults in Santa Barbara, California; carbohydrate entered in the dataset, no text). Step 3: the model was transferred between datasets. 95% confidence intervals (CI) were obtained by cluster BCa bootstrap (2,000 resamples). Results. On the 100 real meals, the carbohydrate mean absolute error (MAE) of the fine-tuned model was 12.2 [9.9, 15.2] g, almost half that of the table (22.6 g; paired difference -10.5 [-16.8, -6.0] g) and nearly a third of that of the same model without fine-tuning (34.3 g), with a bias of -10.3 g; on test dishes it was 2.2 [2.0, 2.5] g. Personalised gradient boosting (M2) reached a correlation between predicted and observed rises of r = 0.63 [0.55, 0.69] in CGMacros (1,286 meals, 45 adults) and 0.49 [0.42, 0.54] in ShanghaiT2DM (2,339 meals, 99 people with T2D), about twice the linear reference (r = 0.32 and 0.23); for a rise above 50 mg/dL, its area under the ROC curve (AUROC) was 0.80 and 0.74, versus 0.64 and 0.68 with pre-meal glucose alone. Adding carbohydrate to pre-meal glucose changed r by 0.095 [0.068, 0.128] in CGMacros and by 0.004 [-0.010, 0.017] in ShanghaiT2DM; difference between datasets: 0.091 [0.060, 0.126]. Applied from one dataset to the other, the model lost correlation ({Delta}r = -0.11 [-0.17, -0.05] towards ShanghaiT2DM and -0.12 [-0.19, -0.04] towards CGMacros). The full "text [->] LLM [->] glucose" chain was evaluated in ShanghaiT2DM only. Conclusion. An open 4-billion-parameter LLM, fine-tuned with QLoRA on public data and usable offline, almost halves the error of carbohydrate counting for real meals described in words compared with an automatically matched table, with an underestimation of large meals. Personalised gradient boosting, based on pre-meal glucose and the person's usual response, doubles the correlation of a linear reference for predicting the glucose rise in both cohorts. Carbohydrate improved this prediction in CGMacros but not in ShanghaiT2DM; the two datasets differ in many respects, and the cause of this gap is not established.
Lombardi, G.; Leung, C. C.-Y.; Burge, R.; Zelaya, F.; Bhattacharyya, S.; Velayudhan, L.
Show abstract
Neuropsychiatric symptoms (NPS) are highly prevalent in Mild Cognitive Impairment (MCI), conceptualized as an intermediate clinical state between normal aging and dementia. Although neuroimaging studies have reported associations between NPS and brain alterations in MCI, longitudinal evidence on NPS and their relationship with multimodal MRI measures remains limited. We aimed to investigate (1) longitudinal changes in NPS, volumetric assessment and measurement of regional cerebral blood flow (CBF) within the medial temporal lobe (MTL); (2) associations between longitudinal changes in MRI measures and NPS; (3) whether baseline MRI measures predicted longitudinal changes in NPS. Thirty individuals with MCI aged over 65 years were recruited, of whom eighteen followed up for a mean period of 29 {+/-} 8.9 months. There were significant reductions in total grey matter and white matter volumes, together with volume loss in the bilateral hippocampus, bilateral parahippocampal gyrus and right entorhinal cortex (pFDR[≤]0.05) over the follow up period. Volume loss in the right parahippocampal gyrus was associated with increase in aggression (pFDR<0.01), and CBF increase in the right hippocampus and parahippocampal gyrus was associated with worsening in agitation (pFDR<0.01). Furthermore, lower baseline volume of the right entorhinal cortex was associated with increase in aggression, and lower baseline CBF in the right hippocampus and parahippocampal gyrus was associated with worsening of agitation (pFDR<0.05). Our findings highlight the importance of identifying NPS in individuals with MCI and the role of structural and perfusion MRI measures as promising biomarkers for monitoring disease progression and improving risk prediction of conversion to dementia.
Ling, Y.; Tsung, S.; Alarid-Escudero, F.; Goldhaber-Fiebert, J. D.; Brandeau, M. L.; Reitsma, M. B.
Show abstract
People experiencing homelessness (PEH) have substantially lower quality-adjusted life expectancy and substantially higher healthcare spending compared to housed populations. California has the largest population of PEH in the United States, and an estimated ten percent of PEH have co-occurring opioid use disorder (OUD). We developed a lifetime, individual-level microsimulation of the 18,928 PEH with OUD in California to evaluate the health and economic effects of expanding access to housing through a Housing First approach, which does not require OUD treatment or abstinence. Increases in permanent housing and shelter access were based on Los Angeles County's Measure A goals, and we accounted for both housing and healthcare costs. Expanding housing access increased mean age at death by 0.34 years, reduced fatal overdoses for the cohort by 39, and produced an incremental cost-effectiveness ratio of $31,083 per quality-adjusted life-year gained, well below conventional willingness-to-pay thresholds. We provide evidence that housing is a high-value health intervention for PEH with OUD in California. Policymakers should support Housing First programs as part of integrated responses to the intertwined crises of homelessness and opioid use.
Hoffman, K. L.; Weil, A. A.; Wetzler, E. A.; Akinola, B. O.; Beauge, A.; Andy, C.; Adams, A.; Frivold, C.; Cox, S. N.; Ehmen, B.; Hatchie, T.; Clark, E. L.; Drummond, M.; Englund, J.; Han, P. D.; Kuntz, J. L.; Lo, N. C.; Lockwood, C.; Starita, L.; Mularski, R. A.; Naleway, A. L.; Carone, M.; Chu, H. Y.
Show abstract
Background: The relationship between recent viral infection and subsequent susceptibility to, or protection against, other viruses remains unclear, particularly across ages and symptom status. Methods: A prospective cohort study enrolled participants aged 6 months to 49 years in Washington and Oregon, USA, 2022-2024. Weekly nasal swabs were tested for SARS-CoV-2, respiratory syncytial virus (RSV), and influenza A/B; additional swabs were collected during symptomatic illness and were also tested for parainfluenza (PIV) 1-4, endemic human coronaviruses (hCOV), adenovirus, human metapneumovirus, human rhinovirus (hRV), and enteroviruses. Proportional hazards regression models were used to characterize associations between incident symptomatic infection (illness) and history of recent illness due to other viruses (within 1-7, 8-14, 15-30, 31-60, 61-90, 91-180 days), adjusted for geographic region, age, sex, daycare/school/work attendance, household size, and number of children in the household. For SARS-CoV-2, influenza, and RSV, these analyses were repeated to study incidence of infection (including asymptomatic cases). Findings: Of 3,560 participants, 2,969 (83%) had [≥]1 illness over 4,299 person-years. Adjusted hazard ratios (aHRs) in the short term ([≤]2 weeks) after an illness were predominantly below 1, indicating that recent viral illness was often associated with lower risk of subsequent illness with other viruses. Across age groups, recent hRV illness was associated with lower risk of hCOV, PIV, and RSV illness (aHRs 0.07-0.27). SARS-CoV-2 illness was associated with a higher short-term risk of adenovirus and PIV illness (aHRs 3.20-3.33). In children aged 6 months-4 years, recent influenza illness was associated with higher short-term risk of SARS-CoV-2 illness. Including asymptomatic infections generally attenuated these associations. Interpretation: Some viral illnesses may alter short-term susceptibility to subsequent viral illnesses. Interactions may differ by age and symptom status. Funding: US Centers for Disease Control and Prevention; Washington State Department of Health.
Zhigulev, A.; Poujade, F.-A.; Gryzunov, N.; Buyan, A.; Li, H.; Bergman, O.; Liu, M.; Pradhananga, S.; Franco-Cereceda, A.; Swoboda, P.; Kulakovskiy, I. V.; Eriksson, P.; Sahlen, P.; Björck, H. M.
Show abstract
BACKGROUND: The chromosome 9p21.3 locus is among the strongest genetic risk loci for cardiovascular disease (CVD), yet the regulatory mechanisms underlying this association remain poorly understood. Previous vascular disease studies focused predominantly on the role of this locus in smooth muscle cell (SMC) dysfunction, whereas its potential role in endothelial cells (ECs) remains less well characterized. Here, we sought to identify candidate regulatory variants at 9p21.3, their distal target genes, and effector transcripts that contribute to CVD-related phenotypes in ECs. METHODS: Promoter capture Hi-C was used to map chromatin interactions at 9p21.3 in primary human ascending aortic SMCs (n=6), supplementing the previously obtained ECs (n=16) dataset. To support the regulatory relevance of these interactions, expression quantitative trait locus analysis was performed in 136 vascular biopsies. Variant effects were prioritized in silico using AlphaGenome. The regulatory activity of the prioritized region was validated using CRISPR activation. RESULTS: Chromatin interaction mapping identified a complex regulatory architecture at the 9p21.3 locus involving 19 putative regulatory elements and 15 genes across ECs and SMCs. MIR31HG showed the most recurrent distal promoter connectivity with the 9p21.3 region. Consistent with the chromatin interaction map, expression of MIR31HG transcripts was associated with genetic variation at the 9p21.3 locus. In ECs, the lncRNA MIR31HG, but not miR-31, was responsive to hemodynamic flow. Depletion of MIR31HG in ECs was associated with endothelial activation. Epigenetic tracks and deep learning predictions highlighted an active enhancer element at 9p21.3, which was subsequently validated by CRISPR activation and increased MIR31HG expression, confirming its capacity to regulate MIR31HG over a large genomic distance of approximately 538 kb. CONCLUSIONS: Altogether, these data implicate lncRNA MIR31HG as a likely effector of the 9p21.3 CVD locus, linking the cis-regulatory element and heritable genetic variation in the respective region to endothelial activation through long-range enhancer regulation.
Kagawa, E.; Kato, M.; Hidaka, T.; Kunita, E.; Hironobe, N.; Matsui, S.; Kobayashi, Y.; Uotani, Y.; Teramoto, T.; Yoshii, K.; Okumoto, T.; Tani, M.; Dote, K.
Show abstract
Background: Artificial intelligence (AI) may improve electrocardiographic differentiation of Takotsubo syndrome from acute anterior ST-segment elevation myocardial infarction (STEMI), but its impact on human diagnostic performance, the added value of visual explanations, and differences by reader background and experience remain uncertain. Methods: Using 12-lead electrocardiogram (ECG) images from two hospitals, we developed and externally validated an ensemble of three deep learning models after self-supervised pretraining and fine-tuning. The diagnostic task used ECGs recorded within 14 days after symptom onset, regardless of reperfusion status. We conducted a four-phase reader study. Phases 1-3 compared unaided interpretation, a reference sheet, and additional AI assistance. In Phase 4, readers evaluated a new ECG set after randomization to reference-sheet support alone, with AI predictions, or with AI predictions plus gradient-weighted class activation mapping (Grad-CAM) heatmaps. Diagnostic correctness was analyzed using generalized linear mixed-effects models. Exploratory analyses examined differences by cardiologist status and experience. Results: Development and external-validation cohorts, including healthy controls, comprised 832 and 95 participants, respectively. For Takotsubo-STEMI discrimination, external-validation accuracy was 0.914, with areas under the receiver operating characteristic and precision-recall curves of 0.970 and 0.984. Seventy healthcare providers completed the reader study. The estimated probability of correct diagnosis increased from 0.632 without assistance to 0.864 with AI assistance (odds ratio [OR], 3.71; P<0.001). In randomized Phase 4, AI assistance improved accuracy vs reference-sheet support alone (OR, 4.28; 95% CI, 3.08-5.95; P<0.001), and Grad-CAM provided additional benefit (OR, 1.53; 95% CI, 1.04-2.25; P=0.033). Exploratory analyses showed greater improvement among less experienced readers; with AI assistance, non-cardiologists outperformed cardiologists. AI-assisted readers remained less accurate overall than the model alone. Conclusions: An ECG image-based ensemble differentiated Takotsubo syndrome from anterior STEMI and improved reader accuracy. Visual explanations provided additional benefit, while experience-related differences highlight the importance of effective human-AI interaction.
Jakob, L.; Mahayri, T. M.; Koutrouli, I. R. A.; Kekrtova, K.; Novakova, D.; Totkova, T.; Spaniel, F.
Show abstract
Schizophrenia spectrum disorders (SSD) have been linked to immune dysregulation and altered intestinal barrier function, although findings across studies remain inconsistent. The present study investigated circulating markers of intestinal permeability and immune activation in individuals with early-stage SSD and healthy controls and examined their associations with clinical outcomes. Serum concentrations of intestinal fatty acid-binding protein (iFABP), zonulin, soluble CD14 (sCD14), bactericidal/permeability-increasing protein (BPI), and lipopolysaccharide-binding protein (LBP) were measured in 36 patients with SSD at baseline and one-year follow-up and in 36 matched healthy controls. Linear mixed-effects models adjusted for age, sex, and body mass index (BMI) were used to evaluate group differences and longitudinal changes. No biomarker differences between patients and controls, and no longitudinal changes within patients, remained statistically significant after correction for multiple testing. Similarly, biomarker concentrations were not associated with symptom severity, global functioning, or neutrophil-to-lymphocyte ratio (NLR). Exploratory analyses identified a robust positive association between BMI and zonulin concentrations and additional associations among iFABP, LBP, and sCD14. Exploratory mixed-effects models identified a negative non-significant association between iFABP and sCD14 after adjustment for age, sex, and BMI. Overall, the findings do not support a generalized alteration of gut permeability or immune activation markers in SSD. Instead, they suggest that variation in gut-related biomarkers may be more strongly influenced by metabolic and individual factors than by diagnosis or current symptom severity.
Arojo, O.; Jones, A. M.
Show abstract
Background: Epidermal Growth Factor Receptor (EGFR) tyrosine kinase inhibitors (TKIs) are widely used for EGFR-mutated non-small cell lung cancer (NSCLC) but have distinct toxicity profiles across generations. Aims: To investigate relationships between EGFR-TKI physicochemical, pharmacokinetic and pharmacological properties and drug-specific adverse drug reactions (ADRs) using national and global pharmacovigilance data. Methods: ADRs for erlotinib, afatinib and osimertinib were extracted from the UK MHRA Yellow Card scheme and WHO VigiAccess. Disproportionality analysis used reporting odds ratios (RORs). Signals were defined as ROR >2, lower 95% confidence interval >1 and >3 cases. Drug properties were curated from ChEMBL, emc, SmPC and EMA sources. Results: VigiAccess contained 88,662 individual case safety reports (163,619 ADRs) and Yellow Card contained 752 reports (1,814 ADRs). Four System Organ Class (SOC) signals were identified in VigiAccess and three in Yellow Card, with five additional Yellow Card Preferred Term signals. Eye disorders showed the strongest consistent association with erlotinib (VigiAccess ROR 4.33; Yellow Card ROR 4.18). Osimertinib was associated with pneumonitis (ROR 7.67) and pulmonary embolism (ROR 3.94). Afatinib showed psychiatric disorder signals in both databases (ROR 3.2-3.3) and diarrhoea (ROR 2.79). Neoplasm-related signals likely reflected confounding by indication. Conclusions: EGFR-TKIs demonstrate generation-specific pharmacovigilance profiles. Integration with drug properties supports mechanistic plausibility for several signals, while the atypical afatinib psychiatric association warrants further investigation. Findings are hypothesis-generating and require confirmation in controlled observational or mechanistic studies.
Taylor, R. W.; Smith, S. J.; Ait-El-Mkadem Saadi, S.; Leighton, J. K.; Taylor, L. S.; Frey, A. M.; Nelson, M. G.; Castaneda-Tamez, P.; Genin, E. C.; Caswell, R.; Dannoura, A.; Shamsheldin, H. E.; Fragaki, K.; Rouzier, C.; He, L.; Hemingbrough, C. V. Y.; Fassad, M. R.; Fouilhoux, A.; Shuaib, T.; Alharby, E.; Rehsi, P.; Rodenburg, R.; Rinne, T.; Janssen, M. C. H.; Hashem, M. O.; Hamid, H.; Bakur, K.; Alrabee, H.; Mohamed, S.; Alansari, R. M.; Al-Owain, M.; Blakely, E. L.; Winter, M.; Almontashiri, N. A. M.; Morris, A. A.; AlAbdi, L.; Alkuraya, F. S.; Robinson, H. K.; McFarland, R.; Alston, C. L.
Show abstract
Mitochondrial disease arises from variants in either the mitochondrial or nuclear genome that impair mitochondrial oxidative phosphorylation. Complex III (ubiquinol:cytochrome c oxidoreductase, CIII2) is a key enzyme of the mitochondrial respiratory chain, catalyzing electron transfer from reduced coenzyme Q10 (ubiquinol) to complex IV, via cytochrome c. CIII2 has one subunit encoded by mitochondrial DNA, the remainder being nuclear encoded. Individuals with CIII2 deficiency present with significant clinical heterogeneity, on account of pathogenic variants in genes encoding structural CIII2 subunits or assembly factors. UQCC1 (GenBank: NM_018244.5) encodes the mitochondrial complex III assembly factor Ubiquinol-cytochrome c reductase complex assembly factor 1, which is currently not associated with human disease. Here we describe nine subjects from eight unrelated families harboring biallelic UQCC1 variants with heterogeneous, early-onset mitochondrial disease presentations and lactic acidosis. Subject-derived fibroblasts demonstrated a decrease in the steady-state levels of UQCC1 and its interacting partner UQCC2, with structural modelling supportive of variant pathogenicity. BN-PAGE and quantitative proteomics demonstrated decreased CIII2 abundance whilst mass spectrometry (MS)-complexome profiling corroborated these findings, identifying impaired respiratory supercomplex formation and a delay in the incorporation of newly synthesized CIII2 subunits into native assemblies. Lentiviral transduction experiments in subject-derived fibroblasts restored the levels of both UQCC1 protein and CIII2 activity. Collectively, these findings provide convergent genetic, biochemical, proteomic and functional evidence that bi-allelic UQCC1 variants disrupt CIII2 assembly, establishing UQCC1 deficiency as a cause of mitochondrial CIII2 deficiency.
Nitola-Pulido, L.; Bonilla, D. S.; Herrera Murcia, E.; Reyes Gavilan, P. A.; Santamaria-Garcia, H.
Show abstract
Exposure to armed conflict can leave persistent effects on mental health, yet its neurophysiological correlates remain poorly characterized, particularly in populations exposed to prolonged violence outside narrowly defined post traumatic stress disorder samples. Here we examined resting state electroencephalographic activity in 122 Colombian adults, including 74 participants exposed to armed conflict and 48 controls. Eyes closed EEG was analysed for 240 seconds per participant. We quantified absolute and relative spectral power across delta, theta, alpha and beta bands and regional functional connectivity using weighted phase lag index in theta, alpha and beta. Group comparisons were adjusted for age, education, sex and assessment year, with correction across 78 EEG contrasts. Median occipital and global alpha power was lower in exposed participants, although adjusted estimates were close to the null and confidence intervals were broad. Temporal theta relative power was lower in exposed participants at nominal significance, but this association did not survive Holm correction. No spectral power or regional connectivity contrast remained significant after multiple testing correction. Exploratory associations between EEG measures and MoCA performance were modest in participants with confirmed scores. These findings provide a systematic characterization of resting EEG after exposure to armed conflict in a population shaped by prolonged violence and displacement. They also indicate that large or widespread alterations in resting spectral power and connectivity were not detectable with the current sample and analytical framework, while defining candidate signals and methodological priorities for larger studies of conflict related brain health.
Falcon-Gracia, A.; Vespalcova, H.; Peligri, D.; Aguilar-Lacasana, S.; Pantano, L.; Nurtdinov, R. N.; Garcia, Z.; Dominguez, A.; Gomez-Herrera, L.; Sanchez, O.; Gascon, M.; Foraster, M.; Guigo, R.; Llurba, E.; Gomez-Roig, M. D.; Rivas, I.; Dadvand, P.; Sunyer, J.; Bustamante, M.; Cosin-Tomas, M.
Show abstract
Fetal development is highly sensitive to environmental influences, including maternal smoking during pregnancy (MSDP) and air pollution. The placenta may play a central role in mediating these effects through epigenetic mechanisms, including small non-coding RNAs (sncRNAs). While most studies have focused on microRNAs (miRNAs), other sncRNA classes remain poorly characterized. Here, we characterized the placental sncRNA landscape and investigated its associations with prenatal exposure to MSDP and air pollution. The sncRNA expression profile of 390 chorionic villi samples from the Barcelona Life Study Cohort (BiSC) was measured using small RNA sequencing and the seqBuster/seqCluster pipeline. MSDP was assessed via questionnaire, and prenatal exposure to black carbon (BC), particulate matter [≤]2.5 m (PM2.5), and nitrogen dioxide (NO2) was estimated using hybrid dispersion-land-use regression (LUR) spatiotemporal models. Associations were evaluated using linear regressions adjusted for relevant covariates. We identified 979 placental sncRNA clusters, assigned primarily to miRNAs (33.2%), piRNAs (26.4%), snoRNAs (12.2%), tRFs (7.8%), and others (20.4%). MSDP was associated with the downregulation of miR-421 and the upregulation of multiple snoRNAs from the imprinted DLK1-DIO3 locus (SNORD113-4, SNORD114-2, SNORD114-3, and SNORD114-17), as well as the piR-26394. No sncRNAs remained significantly associated with air pollution after multiple-testing correction, although several miRNAs showed suggestive associations. Overall, this study provides a comprehensive characterization of placental sncRNAs and identifies previously unreported smoking-associated sncRNAs within the imprinted DLK1-DIO3 locus, expanding our understanding of placental epigenetic responses to prenatal environmental exposures.
Dörfel, R. P.; Ozenne, B.
Show abstract
Estimating the correlation between subject-specific rates of change is a common objective in longitudinal studies. A natural approach is to fit ordinary least squares regressions per subject and correlate the resulting slope estimates directly. Under complete data, this slopes-first estimator is algebraically equivalent to a contrast-based plug-in estimator derived from the joint variance-covariance matrix of the two signals across time. Under independent dropout, however, the slopes-first approach implicitly applies heterogeneous weights across missing-data patterns, introducing bias whenever the cross-covariance between the two signals is not constant over time. We formalise this result and propose an unbiased alternative that estimates the joint covariance matrix before applying the contrast. We further extend our results to include covariates as well. Simulations confirm that the slopes-first estimator is biased under non-constant cross-covariance while the proposed alternative remains unbiased. Further, we show in a real-world example that the two estimators can lead to different scientific conclusions and recommend that researchers working with longitudinal data inspect the missing data pattern to make an informed decision about the estimator of choice.
Dell'Orco, A.; Göschel, L.; Riemann, L. T.; Aydin, S.; Flöel, A.; Tietze, A.; Scheel, M.; Fillmer, A.
Show abstract
Introduction: This study aimed to characterize the neurochemical profile of Alzheimer s disease (AD) continuum stages using proton magnetic resonance spectroscopy (1H-MRS) data for brain metabolites and macromolecules. We used Orthogonal Projection to Latent Structures Discriminant Analysis (OPLS-DA) to classify diagnostic groups and assess feature importance. Methods: Our dataset consisted of 7 T short-echo time 1H-MRS spectra of the posterior cingulate cortex/precuneus of individuals diagnosed with AD (n = 27) or Mild Cognitive Impairment (MCI, n = 27) and Healthy Controls (HC, n = 31). Classification used 100-round Monte Carlo cross-validation based on metabolite concentrations and macromolecular peak intensities. Variable importance in projection (VIP) was assessed across classification models to identify the most relevant neurochemical features. The resulting metabolite rankings were compared with previous literature to evaluate their biological plausibility. Results: Classification was strongest for AD versus HC (AUC=0.88) and lower for comparisons involving MCI. Metabolites contributing most strongly to classification included myo-inositol, ascorbate, and choline-containing compounds, consistent with previous studies on AD-related neurochemical alterations. In addition, specific macromolecular components were identified as relevant features, suggesting a potential contribution of MM signals to the characterization of AD-related neurochemistry. Discussion: These findings support the use of short-echo-time 1H-MRS as a complementary tool for characterizing neurochemical changes in the AD continuum, with potential applications in disease characterization and monitoring, warranting validation in larger cohorts.